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Abstract — We  formulate  and  study  the  thinnest  path  problem 
for  secure  communication  in  wireless  ad  hoc  networks.  The  objec¬ 
tive  is  to  find  a  path  from  a  source  to  its  destination  that  results  in 
the  minimum  number  of  nodes  overhearing  the  message  by  a  judi¬ 
cious  choice  of  relaying  nodes  and  their  corresponding  transmission 
powers.  We  adopt  a  directed  hypergraph  model  of  the  problem  and 
establish  the  NP-completeness  of  the  problem  in  2-D  networks.  We 
then  develop  two  polynomial-time  approximation  algorithms  that 
offer  and  =  approximation  ratios  for  general  directed 
hypergraphs  (which  can  model  nonisotropic  signal  propagation  in 
space)  and  constant  approximation  ratios  for  ring  hypergraphs 
(which  result  from  isotropic  signal  propagation).  We  also  consider 
the  thinnest  path  problem  in  1-D  networks  and  1-D  networks  em¬ 
bedded  in  a  2-D  field  of  eavesdroppers  with  arbitrary  unknown  loca¬ 
tions  (the  so-called  1 .5-D  networks).  We  propose  a  linear-complexity 
algorithm  based  on  nested  backward  induction  that  obtains  the  op¬ 
timal  solution  for  both  1-D  and  1.5-D  networks.  This  algorithm  does 
not  require  the  knowledge  of  eavesdropper  locations  and  achieves 
the  best  performance  offered  by  any  algorithm  that  assumes  com¬ 
plete  location  information  of  the  eavesdroppers. 

Index  Terms — Approximation  algorithms,  approximation  ratio, 
hypergraph,  NP-complete,  secure  communication,  thinnest  path. 


I.  Introduction 
A.  Thinnest  Path  Problem 

IN  THIS  paper,  we  consider  the  thinnest  path  problem  for 
secure  communication  in  wireless  ad  hoc  networks.  For  a 
given  source  and  a  destination,  the  thinnest  path  problem  asks 
for  a  path  from  the  source  to  the  destination  that  results  in  the 
minimum  number  of  nodes  hearing  the  message.  Such  a  path 
is  achieved  by  carefully  choosing  a  sequence  of  relaying  nodes 
and  their  corresponding  transmission  powers. 

At  first  glance,  one  may  wonder  whether  the  thinnest  path 
problem  is  simply  a  shortest  path  problem  with  the  weight  of 
each  hop  given  by  the  number  of  nodes  that  hear  the  message 
in  that  hop.  Realizing  that  a  node  may  be  within  transmission 
range  of  multiple  relaying  nodes  and  should  not  be  counted  mul¬ 
tiple  times  in  the  total  weight  (referred  to  as  the  width)  of  the  re¬ 
sulting  path,  we  see  that  the  thinnest  path  problem  does  not  have 
a  simple  cost  function  that  is  summable  over  edges.  Rather,  the 
width  of  a  path  is  given  by  the  cardinality  of  the  union  of  all 
receiving  nodes  in  each  hop,  which  is  a  highly  nonlinear  func¬ 
tion  of  the  weight  of  each  hop.  One  may  then  wonder  whether 

Manuscript  received  October  22,  2013;  revised  February  28,  2014;  accepted 
April  12,  2014;  approved  by  IEEE/ACM  Transactions  on  Networking  Ed¬ 
itor  G.  Xue.  Date  of  publication  June  03,  2014;  date  of  current  version  August 
14,  2015.  This  work  was  supported  by  the  Army  Research  Laboratory  Network 
Science  CTA  under  Cooperative  Agreement  W911NF-09-2-0053. 

J.  Gao  and  Q.  Zhao  are  with  the  University  of  California,  Davis,  Davis,  CA 
95616  USA  (e-mail:  jhgao@ucdavis.edu;  qzhao@ucdavis.edu). 

A.  Swami  is  with  the  Army  Research  Laboratory,  Adelphi,  MD  20783  USA 
(e-mail:  a.swami@ieee.org). 

Color  versions  of  one  or  more  of  the  figures  in  this  paper  are  available  online 
at  http://ieeexplore.ieee.org. 

Digital  Object  Identifier  10.1109/TNET.2014.2321159 


we  can  redefine  the  weight  of  each  hop  as  the  number  of  nodes 
that  hear  the  message  for  the  first  time.  Such  a  definition  of  edge 
weight  indeed  leads  to  a  summable  cost  function.  Unfortunately, 
in  this  case,  the  edge  weight  cannot  be  predetermined  until  the 
thinnest  path  from  the  source  to  the  destination  in  question  has 
already  been  established. 

A  more  fundamental  difference  between  the  thinnest  path  and 
the  shortest  path  problems  is  that  the  thinnest  paths  from  a  single 
source  to  all  other  nodes  in  the  network  do  not  form  a  tree.  In 
other  words,  the  thinnest  path  to  a  node  does  not  necessarily  go 
through  the  thinnest  path  to  any  of  its  neighbors.  The  loss  of  the 
tree  structure  is  one  of  the  main  reasons  that  the  thinnest  path 
problem  is  much  more  complex  than  the  shortest  path  problem. 
Indeed,  as  shown  in  this  paper,  the  thinnest  path  problem  is 
NP-complete,  which  is  in  sharp  contrast  with  the  polynomial 
nature  of  the  shortest  path  problem. 

Another  aspect  that  complicates  the  problem  is  the  choice 
of  the  transmission  power  at  each  node  (within  a  maximum 
value  that  may  vary  across  nodes).  In  this  case,  the  network 
cannot  be  modeled  as  a  simple  graph  in  which  the  neighbors 
of  each  node  are  prefixed.  In  this  paper,  we  adopt  the  directed 
hypergraph  model  that  easily  captures  the  choice  of  different 
neighbor  sets  (corresponding  to  different  transmission  powers) 
at  each  node.  While  a  graph  is  given  by  a  vertex  set  V  and 
an  edge  set  E  consisting  of  cardinality-2  subsets  of  V,  a  hy¬ 
pergraph  [1]  is  free  of  the  constraint  on  the  cardinality  of  an 
edge.  Specifically,  any  nonempty  subset  of  V  can  be  an  element 
(referred  to  as  a  hyperedge)  of  the  edge  set  E.  Hypergraphs 
can  thus  capture  group  behaviors  and  higher-dimensional  re¬ 
lationships  in  complex  networks  that  are  more  than  a  simple 
union  of  pairwise  relationships  [2].  In  a  directed  hypergraph  [3], 
each  hyperedge  is  directed,  going  from  a  source  vertex  to  a 
nonempty  set  of  destination  vertices.  An  example  is  given  in 
Fig.  1(a)  where  we  have  two  directed  hyperedges  rooted  at  a 
source  node  v  with  each  hyperedge  modeling  a  neighbor  set 
of  v  under  a  specific  power.  The  directed  hypergraph  model  of 
the  thinnest  path  problem  is  thus  readily  seen:  Rooted  at  each 
node  are  multiple  directed  hyperedges,  each  corresponding  to 
a  distinct  neighbor  set  feasible  under  the  maximum  transmis¬ 
sion  power  of  this  node.  The  problem  is  then  to  find  a  min¬ 
imum-width  hyperpath  from  the  source  to  the  destination  where 
the  width  of  a  hyperpath  is  given  by  the  cardinality  of  the  union 
of  the  hyperedges  on  this  hyperpath. 

B.  Main  Results 

Based  on  the  directed  hypergraph  formulation,  we  show  that 
the  thinnest  path  problem  in  2-D  networks  is  NP-complete 
even  under  a  simple  disk  propagation  model.  This  result  is 
established  through  a  reduction  from  the  minimum  dominating 
set  (MDS)  problem  in  graphs,  a  classic  NP-complete  problem. 
The  most  challenging  part  of  this  reduction  is  to  show  the 
reduced  problem  is  realizable  under  a  2-D  disk  model  that 
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Fig.  1 .  (a)  Directed  hypergraph  (two  hyperedges  represented  by  the  green  solid 
line  and  red  dashed  lines,  respectively),  (b)  Disk  hypergraph  (three  hyperedges 
represented  by  the  green  solid  line,  red  solid  line,  and  red  dashed  lines,  respec¬ 
tively).  (c)  Unit  disk  hypergraph  (four  hyperedges  represented  by  the  green  solid 
line,  blue  solid  line,  red  solid  line,  and  red  dashed  lines,  receptively). 


has  specific  geometrical  properties  that  need  to  be  preserved 
in  the  reduction.  We  further  establish  that  even  with  a  fixed 
transmission  power  at  each  node  (in  this  case,  the  resulting 
hypergraph  degenerates  to  a  standard  graph),  the  thinnest  path 
problem  is  NP-complete.  We  then  propose  two  polynomial-time 
approximation  algorithms  that  offer  and  =  approx¬ 
imation  ratios  for  general  directed  hypergraphs  (which  can 
model  nonisotropic  signal  propagation  in  space)  and  constant 
approximation  ratios  for  ring  hypergraphs  (which  result  from 
isotropic  signal  propagation).  Here,  n  is  the  total  number  of 
vertices. 

We  also  establish  the  polynomial  nature  of  the  problem  in 
1-D  and  1.5-D  networks,  where  a  1.5-D  network  is  a  1-D  net¬ 
work  embedded  in  a  2-D  field  of  eavesdroppers  with  arbitrary 
unknown  locations.  We  propose  an  algorithm  based  on  a  nested 
backward  induction  (NBI)  starting  at  the  destination.  We  show 
that  this  NBI  algorithm  has  0(n)  time  complexity.  Since  the 
size  of  the  input  data  is  0(n),  the  proposed  algorithm  is  order- 
optimal.  It  solves  the  thinnest  path  problem  in  both  the  1-D  and 
1.5-D  networks.  In  particular,  no  algorithm,  even  with  complete 
location  information  of  the  eavesdroppers,  can  obtain  a  thinner 
path  than  the  NBI  algorithm,  which  does  not  require  knowledge 
of  eavesdropper  locations. 

In  a  broader  context,  the  concepts  and  techniques  of  directed 
line  crossing  and  exposed  disk  hypergraphs  introduced  in  this 
paper  for  preserving  geometrical  properties  when  establishing 
the  NP-completeness  of  the  problem  provide  new  tools  for 
complexity  studies  in  geometrical  hypergraphs  and  graphs.  The 
bounding  techniques  and  the  use  of  sphere  packing  results  in 
analyzing  the  performance  of  the  two  approximation  algorithms 
may  also  find  other  applications  in  algorithmic  analysis. 

In  the  context  of  secure  communications,  the  motivation  for 
the  thinnest  path  problem  is  to  reduce  the  risk  of  information 
leakage  by  minimizing  the  number  of  in-network  nodes  and  as 
well  as  eavesdroppers  overhearing  messages  that  are  intended 
only  for  a  specific  destination  node.  The  problem  may  also  have 
implications  from  the  energy  efficiency  perspective.  Nodes  that 
receive  a  signal  may  attempt  to  decode  it,  even  if  they  are  not 
in  the  optimal  relay  path.  This  may  be  particularly  important 
in  a  duty-cycled  sensor  network  where  inadvertent  signals  may 
wake  up  sensors  and  cause  unnecessary  energy  consumption. 

C.  Related  Work 

There  is  a  large  body  of  literature  on  security  issues  in  wire¬ 
less  ad  hoc  networks  (see,  for  example,  [4]  and  [5]).  However, 
the  thinnest  path  problem  has  not  been  studied  in  the  literature 


except  in  [6].  Chechik  et  al  studied  the  thinnest  path  (referred 
to  as  the  secluded  path  in  [6])  and  the  thinnest  Steiner  tree  in 
graphs.  They  showed  that  the  problem  in  a  general  graph  is 
NP-complete  and  strongly  inapproximable.  They  proposed  an 
algorithm  with  an  approximation  ratio  of  \/A  +  3  for  bounded- 
degree  graphs  where  A  is  the  maximum  degree.  They  further 
studied  the  problem  in  several  special  graph  models  including 
hereditary  graphs  and  planar  graphs.  However,  their  study  fo¬ 
cuses  on  the  problem  in  topological  graphs,  whereas  we  focus  on 
hypergraphs  and  geometric  graphs.  The  complexity  results  ob¬ 
tained  in  [6]  do  not  apply  to  special  hypergraphs  satisfying  cer¬ 
tain  geometric  properties  that  result  naturally  from  the  commu¬ 
nication  problem  studied  in  this  paper.  This  paper  also  includes 
several  new  complexity  results  on  the  thinnest  path  problem 
under  the  geometrical  graph  models.  Specifically,  we  establish 
the  NP-completeness  of  the  problem  in  2-D  disk  graphs  and 
3-D  unit  disk  graphs.  Furthermore,  [6]  and  this  paper  use  dif¬ 
ferent  techniques  in  the  complexity  analysis.  In  particular,  [6] 
demonstrates  the  NP-completeness  by  constructing  a  reduction 
from  the  red-blue  vertex  cover  problem,  whereas  we  construct 
a  reduction  from  the  MDS  problem.  The  reason  that  different 
techniques  are  needed  is  that  [6]  focuses  on  topological  graphs 
and  topological  relationships  (i.e.,  who  is  connected  to  whom) 
that  are  easier  to  maintain  during  the  reduction.  In  our  case,  we 
consider  geometrical  models  that  dictate  not  only  who  is  con¬ 
nected  to  whom,  but  also  the  relative  positions  (e.g.,  connecting 
edges  cannot  arbitrarily  cross  each  other  without  consequences). 
In  order  to  preserve  all  the  geometrical  properties  of  the  orig¬ 
inal  network,  the  reduction  needs  to  be  carefully  constructed.  In 
particular,  our  techniques  of  using  directed  line  crossing  and  ex¬ 
posed  disk  hypergraphs  are  novel  concepts  for  maintaining  ge¬ 
ometrical  properties.  The  results  in  [6]  and  this  work  thus  com¬ 
plement  each  other  to  provide  a  more  complete  picture  of  the 
thinnest  path  problem  under  different  (hyper)graph  models. 

The  shortest  path  problem  in  hypergraphs  remains  a  polyno¬ 
mial-time  problem  as  its  counterpart  under  the  graph  model. 
Existing  work  on  both  the  static  and  dynamic  version  of  the 
shortest  path  problem  in  hypergraphs  can  be  found  in  [3]  and 
[7]-[9].  As  discussed  earlier,  the  thinnest  path  problem  is  fun¬ 
damentally  different  and  significantly  more  complex  than  the 
shortest  path  problem. 

The  widest  path  problem  has  been  well  studied  under  the 
graph  model  [10],  [11],  and  these  existing  results  can  be  easily 
extended  to  hypergraphs.  The  widest  path  problem  asks  for  a 
path  whose  minimum  edge  weight  along  the  path  is  maximized. 
In  other  words,  the  width  of  a  path  is  given  by  the  minimum 
edge  weight  on  that  path,  which  is  different  from  the  definition 
of  path  width  in  the  thinnest  path  problem  studied  in  this  paper. 
As  a  consequence,  the  widest  path  problem  is  not  the  comple¬ 
ment  of  the  thinnest  path  problem.  Since  the  tree  structure  is 
preserved  in  the  widest  path  problem  (i.e.,  the  widest  path  to  a 
node  must  go  through  the  widest  path  to  one  of  its  neighbors),  it 
remains  a  polynomial-time  problem.  The  thinnest  path  problem, 
however,  is  NP-complete  in  general. 

Another  related  problem  is  topology  control,  where  the  ob¬ 
jective  is  to  design  the  transmission  power  of  each  node  such 
that  the  maximum  interference  in  the  network  (measured  by  the 
maximum  in-degree  over  all  nodes)  is  minimized  under  the  con¬ 
straint  that  the  network  is  connected.  Rickenbach  et  al.  [12] 
and  Halldorsson  et  al.  [13]  studied  the  problem  in  1-D  and 
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2-D  wireless  ad  hoc  networks.  The  focus  of  [12]  and  [13]  is 
on  developing  approximation  algorithms;  the  hardness  of  the 
problem  remains  open.  While  both  the  thinnest  path  and  the 
topology  control  problems  involve  the  design  of  transmission 
powers,  the  objectives  are  fundamentally  different.  As  a  result, 
they  call  for  different  techniques  in  both  complexity  analysis 
and  algorithm  design. 

In  the  general  context  of  algorithmic  studies  in  hypergraphs, 
Ausiello  et  al.  [14]  tackled  the  problem  of  finding  the  fi- op¬ 
timal  hyperpath  where  fi  is  a  general  measure  on  hyperpaths 
that  satisfies  a  certain  monotone  property.  They  established 
the  NP-completeness  of  this  problem  for  general  measures. 
The  thinnest  path  problem  can  be  seen  as  a  //-optimal  traversal 
problem  with  the  measure  \i  given  by  the  number  of  vertices 
covered  by  the  path.  Since  this  is  a  special  measure,  their  NP-com¬ 
pleteness  result  developed  under  general  measures  does  not 
apply.  Furthermore,  in  many  applications,  the  resulting  hyper¬ 
graphs  have  certain  topological  and/or  geometrical  properties, 
and  the  computational  complexities  under  these  special  models 
require  separate  analysis. 

II.  Problem  Formulation 

A.  Basic  Concepts  of  Directed  Hypergraphs 

A  directed  hypergraph  H  =  (V,  E)  consists  of  a  set  V  of 
vertices  and  a  set  E  of  directed  hyperedges  [3].1  Each  directed 
hyperedge  e  £  E  has  a  single  source  vertex  se  and  a  nonempty 
set  of  destination  vertices  Te .  We  let  n  =  \V\  denote  the  number 
of  vertices. 

A  disk  hypergraph  is  a  special  directed  hypergraph  whose 
topology  is  determined  by  a  set  of  points  V  =  {ui, . . . ,  vn} 
located  in  a  d-dimensional  Euclidean  space  and  a  maximum 
range  Ri  associated  with  each  vertex  V{ .  There  exists  a  hyper¬ 
edge  e  from  source  se  =  Vi  to  destination  set  Te  if  and  only  if 
Te  consists  of  vertices  located  within  the  d-dimensional  sphere 
centered  at  Vi  with  a  radius  r  £  (0,  i2*].  A  unit  disk  hyper¬ 
graph  (UDH)  is  a  disk  hypergraph  with  unit  maximum  range 
(Ri  =  1)  for  all  vertices.  Fig.  1  shows  examples  of  a  directed 
hypergraph,  a  disk  hypergraph,  and  a  unit  disk  hypergraph. 

A  ring  hypergraph  is  a  generalized  disk  hypergraph  where 
associated  with  each  vertex  V{  is  a  minimum  range  as  well 
as  a  maximum  range  R{.  Hyperedges  rooted  at  V{  are  formed 
by  spheres  centered  at  V{  with  radii  satisfying  V{  <  r  <  Ri .  It 
is  easy  to  see  that  a  disk  hypergraph  is  a  ring  hypergraph  with 
ri  =  0,  a  disk  graph  is  a  ring  hypergraph  with  r*  =  Ri  for 
all  z,  and  a  unit  disk  graph  (UDG)  is  a  ring  hypergraph  with 
rL  =  Rj  =  1  for  all  z. 

B.  Thinnest  Path  Problem 

Consider  a  wireless  ad  hoc  network  with  n  nodes  located 
in  a  d-dimension  Euclidean  space.  Each  node  can  choose  the 
power,  within  a  maximum  value,  for  the  transmission  of  each 
message.  The  chosen  power,  along  with  the  signal  propagation 
model,  determines  the  set  of  neighbors  that  can  hear  the  mes¬ 
sage.  The  maximum  transmission  power  is  in  general  different 
across  nodes.  The  objective  is  to  find  a  path  between  a  given 
source-destination  pair  that  involves  the  minimum  number  of 
nodes  hearing  the  message. 

As  discussed  in  Section  I- A,  we  formulate  the  problem  using 
a  directed  hypergraph.  Each  node  is  a  vertex.  The  directed  hy- 

Tn  [3],  it  is  referred  to  as  the  forward  hyperarcs. 


TABLE  I 
Notations 

e  :  hyperedge  se  :  source  of  e  Te  :  destinations  of  e 
L  :  hyperpath  L  :  cover  of  L  W  ( L )  :  width  of  L 

peredges  rooted  at  a  node  are  given  by  distinct  neighbor  sets 
of  this  node  feasible  under  its  maximum  transmission  power 
and  the  signal  propagation  model.  Under  a  general  nonisotropic 
propagation  model,  we  end  up  with  a  general  hypergraph.  The 
only  property  the  resulting  hypergraph  has  is  the  monotonicity 
of  the  hyperedge  set.  Specifically,  the  hyperedges  rooted  at  each 
node  can  be  ordered  in  such  a  way  (say,  a,  e2, . . . ,  ei)  that 
Tei  C  Tei+ 1  and  \Te.  |  =  \Tei+1 1  -  1.  This  is  due  to  the  nature 
of  wireless  broadcasting  where  nodes  reachable  under  transmis¬ 
sion  power  rj  can  also  be  reached  under  any  power  greater  than 
rj.  Under  an  isotropic  propagation  model,  we  end  up  with  a  disk 
hypergraph.  If  all  nodes  have  the  same  maximum  range,2  we 
have  a  unit  disk  hypergraph.  This  hypergraph  model  also  applies 
to  networks  with  eavesdroppers.  Each  eavesdropper  can  be  seen 
as  a  node  with  zero  transmission  range.  It  is  thus  a  vertex  with 
no  outgoing  hyperedges. 

Given  a  source-destination  pair  (s,  t),  a  hyperpath  from  s  to  t 
is  defined  as  a  sequence  of  hyperedges  L  =  {ei,...,em}  such 
that  sei  £  Tei_1  for  1  <  z  <  m,  sei  =  s,  and  t  £  Tem.  Define 
the  cover  L  of  L  to  be  the  set  of  vertices  in  L,  i.e., 

i  =  uf=1reju{Sei}. 

The  width  W(L )  is  then  given  by 

W{L)t\L\. 

The  thinnest  path  problem  asks  for  a  hyperpath  from  s  to  t 
with  the  minimum  width.  Note  that  choosing  a  hyperedge  e  = 
{se,Te}  simultaneously  chooses  the  relaying  node  se  and  its 
transmission  power  (determined  by  Te). 

Table  I  lists  notations  used  throughout  the  paper. 

III.  NP-Complete  Problems 

In  this  section,  we  show  that  the  thinnest  path  (TP)  problem 
is  NP-complete  in  several  special  geometric  hypergraphs  and 
graphs.  This  implies  the  NP-completeness  of  the  problem  in 
general  directed  hypergraphs. 

A.  TP  in  2-D  Disk  Hypergraphs 

In  this  section,  we  prove  the  NP-completeness  of  the  thinnest 
path  problem  in  2-D  disk  hypergraphs.  While  a  stronger  result 
is  shown  in  Section  III-B,  the  proof  of  this  result  provides  the 
main  building  block  for  the  proof  of  the  next  result. 

The  result  is  established  through  a  reduction  from  the 
MDS  [15]  problem.  The  MDS  problem  asks  for  the  minimum 
subset  of  vertices  in  a  given  graph  such  that  every  vertex  in 
the  graph  is  either  in  the  subset  or  a  direct  neighbor  of  a  vertex 
in  the  subset.  The  following  theorem  formally  establishes  the 
polynomial  reduction  (denoted  by  <p)  from  MDS  to  TP  in  2-D 
disk  hypergraphs.  Since  the  thinnest  path  problem  is  clearly  in 
the  NP  space,  this  theorem  establishes  the  NP-completeness  of 
TP  in  2-D  disk  hypergraphs. 

Theorem  1:  MDS  <p  TP  in  2-D  disk  hypergrpahs. 

To  prove  Theorem  1 ,  consider  an  MDS  problem  in  an  arbitrary 
graph  G.  We  first  construct  a  general  directed  hypergraph  Hi 

transmission  range  and  transmission  power  are  used  interchangeably. 


GAO  et  al. :  THINNEST  PATH  PROBLEM 


1179 


based  on  G  such  that  a  thinnest  path  in  Hi  leads  to  an  MDS  in 
G.  The  main  challenge  in  the  proof  is  to  show  that  Hi  is  realizable 
under  a  2-D  disk  model.  There  are  two  main  difficulties.  First, 
line  crossing  is  inevitable  when  we  draw  Hi  on  a  2-D  plane.  The 
implementation  of  hyperedges  that  cross  each  other  needs  special 
care  to  avoid  unwanted  overhearing  that  may  render  the  reduction 
invalid.  Second,  the  geometric  structure  of  2-D  disk  hypergraphs 
dictates  that  there  are  at  most  five  vertices  (even  with  arbitrary 
ranges)  that  can  reach  a  common  sixth  vertex  but  not  each  other. 
It  is  thus  challenging  to  implement  a  vertex  with  up  to  n  incoming 
hyperedges  in  Hi  while  preserving  the  reduction. 

Our  main  approach  to  overcoming  the  above  difficulties  is  to 
allow  directed  overhearing.  Specifically,  messages  transmitted 
along  one  hyperedge  may  be  heard  by  vertices  implementing  an¬ 
other  hyperedge  in  Hi ,  but  not  vice  versa.  By  carefully  choosing 
the  directions  of  the  introduced  overhearing,  we  ensure  that  the 
resulting  2-D  disk  hypergraph  H2 ,  while  having  a  different  topo¬ 
logical  structure  fromfii ,  preserves  the  reduction  fromMD  S  in  G. 

Another  challenge  in  constructing  H2  is  to  ensure  the  polyno¬ 
mial  nature  of  the  reduction.  The  number  of  additional  vertices 
added  in  H2  needs  to  be  in  a  polynomial  order  in  the  size  of  G. 
This  often  limits  the  use  of  reduced  transmission  ranges  as  a  way 
to  avoid  unwanted  overhearing:  Exponentially  small  transmis¬ 
sion  ranges  may  require  exponentially  many  vertices  to  connect 
two  fixed  points. 

A  detailed  proof  is  given  in  Appendix  A. 

B.  TP  in  2-D  Unit  Disk  Hypergraphs 

We  now  establish  the  NP-completeness  of  TP  in  2-D  unit  disk 
hypergraphs.  The  proof  builds  upon  the  proof  of  Theorem  1 .  The 
only  difference  is  that  when  implementing  the  general  directed 
hypergraph  Hi,  we  no  longer  have  the  freedom  of  choosing  the 
maximum  transmission  range  of  each  vertex.  This  presents  a 
nontrivial  challenge.  As  stated  in  Section  III- A,  our  approach 
to  circumvent  the  constraints  imposed  by  the  geometrical  struc¬ 
tures  of  2-D  disk  hypergraphs  is  to  allow  directed  overhearing, 
which  is  achieved  by  carefully  choosing  different  maximum 
transmission  ranges  of  various  vertices.  To  implement  a  2-D 
UDH  for  the  reduction,  however,  all  vertices  must  have  the  same 
maximum  transmission  range. 

To  address  this  issue,  we  introduce  a  special  type  of  disk  hy¬ 
pergraph,  called  exposed  disk  hypergraphs ,  and  show  that  TP 
in  k- D  exposed  disk  hypergraphs  can  be  reduced  to  TP  in  k- D 
UDH  for  any  k  >  2.  We  then  show  that  the  2-D  disk  hypergraph 
H2  constructed  in  the  proof  of  Theorem  1  can  be  modified  to 
an  exposed  hypergraph  while  preserving  the  reduction.  We  thus 
arrive  at  the  NP-completeness  of  TP  in  2-D  UDH  based  on  the 
transitivity  of  polynomial  time  reduction. 

Definition  1:  In  a  disk  hypergraph  H  =  (V,  E ),  let  rv  denote 
the  closest  nonneighbor3  of  v.  Define4 

€v  —  —  {dipj^Ty)  Ry) 

where  d{v,  tv)  is  the  distance  between  v  and  rv  ( ev  is  set  to  1 
when  v  does  not  have  nonneighbors).  An  exposed  area  $>v  of  v 
is  defined  as  \ 

-^v ,Rv-\-ev  \  |^J 

\uev 

where  Dv,r  denotes  the  closed  ball  centered  at  v  with  radius  r. 
A  disk  hypergraph  is  exposed  if  every  vertex  has  a  nonempty 

3  A  vertex  is  a  nonneighbor  of  v  if  it  is  outside  the  maximum  range  Rv  of  v. 

4The  parameter  |  can  be  changed  to  an  arbitrary  positive  value  smaller  than  1 . 


Fig.  2.  Exposed  hypergraphs  and  exposed  areas  (H1  is  not  exposed  since  v  has 
an  empty  exposed  area;  H2  and  H3  are  exposed). 


exposed  area.  Fig.  2  demonstrates  one  nonexposed  hypergraph 
and  two  exposed  hypergraphs  with  the  exposed  areas. 

Lemma  1:  TP  in  k- D  exposed  disk  hypergraphs  <P  TP  in 
k- D  UDH. 

Proof:  The  basic  idea  is  to  place  super  vertices  at  specific 
locations  in  exposed  areas  to  force  vertices  on  a  thinnest  path 
to  use  transmission  ranges  smaller  than  the  maximum  value. 
The  problem  is  thus  transformed  to  the  case  with  disk  hyper¬ 
graphs  where  vertices  may  have  different  maximum  transmis¬ 
sion  ranges.  A  detailed  proof  is  given  in  Appendix  B.  ■ 

With  Lemma  1  providing  a  bridge  between  disk  and  unit  disk 
hypergraphs,  all  we  need  to  show  is  that  MDS  can  be  reduced 
to  TP  in  2-D  exposed  disk  hypergraphs. 

Lemma  2:  MDS  <P  TP  in  2-D  exposed  disk  hypergrpahs. 

Proof:  See  Appendix  C.  ■ 

Based  on  Lemmas  1  and  2,  we  arrive  at  the  following 
theorem. 

Theorem  2:  MDS  <P  TP  in  2-D  UDH. 

C.  TP  in  2-D  Disk  Graphs  and  3-D  Unit  Disk  Graphs 

In  this  section,  we  consider  the  thinnest  path  problem  in  disk 
graphs  and  UDGs.  Recall  that  disk  and  unit  disk  graphs  are 
special  ring  hypergraphs  with  r,[  =  R{  and  n  —  R{  — 
respectively.  In  other  words,  they  can  be  seen  as  hypergraphs 
where  each  vertex  has  only  one  outgoing  hyperedge  directed  to 
its  prefixed  neighbor  set  (determined  by  its  fixed  transmission 
power).  This  also  shows  that  disk  hypergraphs  and  disk  graphs 
are  not  special  cases  of  each  other.  Given  the  same  set  of  vertices 
and  their  associated  maximum  ranges,  a  disk  hypergraph  has  a 
topology  different  from  that  of  a  disk  graph:  Each  vertex  in  gen¬ 
eral  has  more  than  one  outgoing  hyperedge  due  to  the  freedom 
of  using  smaller  transmission  ranges.  The  same  holds  for  UDH 
and  UDG.  As  a  consequence,  the  complexity  of  TP  in  disk  and 
unit  disk  graphs  cannot  be  inferred  from  Theorems  1  and  2  and 
needs  to  be  studied  separately. 

Theorem  3:  MDS  <P  TP  in  2-D  disk  graphs. 

Proof:  In  the  proof  of  Theorem  1 ,  the  vertices  along  the 
thinnest  path  in  the  constructed  2-D  disk  hypergraph  H2  all  use 
their  maximum  ranges.  Thus,  MDS  in  G  can  be  reduced  to  TP 
in  a  disk  graph  constructed  from  H2  by  including  only  those 
hyperedges  associated  with  the  maximum  range  of  each  vertex. 

■ 

Next,  we  consider  TP  in  UDG.  Unfortunately,  the  ap¬ 
proach  through  exposed  disk  hypergraphs  used  in  showing 
the  NP-completeness  of  TP  in  UDH  does  not  apply  since  it 
hinges  on  vertices  being  able  to  use  any  transmission  range 
smaller  than  a  maximum  value.  The  difficulty,  however,  can  be 
circumvented  for  3-D  UDG  as  shown  in  the  following  theorem. 

Theorem  4:  MDS  in  degree-3  graphs  <P  TP  in  3-D  UDG. 

The  proof  is  similar  to  that  of  Theorem  1  with  two  main  dif¬ 
ferences.  First,  line  crosses  are  implemented  by  using  the  third 
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dimension  to  “go  around,”  rather  than  using  different  transmis¬ 
sion  ranges  (a  luxury  absent  in  UDG)  to  create  directed  crosses. 
Second,  reduction  from  MDS  in  graphs  with  a  maximum  de¬ 
gree  of  3  ensures  that  there  are  at  most  four  incoming  edges 
to  each  super  vertex  in  the  reduced  UDG.  This  makes  the  geo¬ 
metric  constraint  on  the  number  (at  most  11  in  a  3-D  Euclidean 
space)  of  vertices  that  can  reach  a  common  vertex  but  not  each 
other  inconsequential.5  A  detailed  proof  is  given  in  Appendix  D. 

Note  that  using  a  reduction  from  MDS  in  graphs  with  a 
constant  maximum  degree  rather  than  MDS  in  general  graphs 
leads  to  a  weaker  statement.  While  MDSs  in  both  cases  are 
NP-complete,  the  former  is  approximable  with  a  constant  ratio, 
and  the  latter  a  ratio  of  0(log  n).  Theorems  1-3  thus  give  a 
log  n-order  lower  bound  on  the  approximation  ratio  of  those 
problems,  whereas  Theorem  4  provides  a  constant  lower  bound. 

IV.  Polynomial  Complexity  Problems 

In  this  section,  we  consider  the  thinnest  path  problem  in  1  -D 
networks.  We  show  that  the  problem  is  polynomial  time  by  con¬ 
structing  an  algorithm  with  time  complexity  of  0{n).  Since  the 
input  data  has  size  0(n),  the  proposed  algorithm  is  order-op¬ 
timal.  We  then  consider  the  1.5-D  problem  and  show  that  the 
algorithm  developed  for  1-D  networks  directly  applies  to  the 
1.5-D  problem. 

A.  1-D  Networks 

Consider  a  network  under  a  general  propagation  model  with 
n  nodes  located  on  a  straight  line.  Each  vertex  V{  is  associated 
with  a  coordinate  X{  on  the  line  (the  vertex  index  Vi  and  its  lo¬ 
cation  Xi  are  often  used  interchangeably).  Without  loss  of  gen¬ 
erality,  we  assume  that  x\  <  x2  <  . . .  <  xn. 

It  is  clear  that  every  node  located  between  the  source  s  and  the 
destination  t  (see  Fig.  3)  will  hear  the  message  no  matter  which 
path  is  chosen,  and  all  nodes  to  the  right  of  t  can  be  excluded 
from  the  thinnest  path.  Therefore,  finding  the  thinnest  path  is 
equivalent  to  minimizing  the  number  of  vertices  to  the  left  of  s 
that  can  overhear  the  message.  The  problem  is  nontrivial.  Due 
to  the  arbitrariness  of  the  node  locations  and  propagation  range, 
a  forward  path  (i.e.,  every  hop  moves  the  message  to  the  right 
toward  t)  from  s  to  t  may  not  exist,  and  nodes  to  the  left  of  s 
may  need  to  act  as  relays.  The  question  is  thus  how  to  efficiently 
find  out  whether  a  forward  path  exists  and,  if  not,  which  set  of 
nodes  to  the  left  of  s  need  to  relay  the  message. 

We  propose  an  algorithm  based  on  NBI.  For  each  vertex  v, 
we  define  its  predecessor  pv  to  be  the  nearest  vertex  on  the  left 
side  of  v  that  can  reach  v 

Pv  =  argmax{x„  :  xu  <  xv, 

u£V 

Be  E  E  s.t.  se  =  u  and  v  £  Te}.  (1) 
Thus,  in  order  to  reach  v,  its  predecessor  pv  or  a  vertex  to  the 
left  of  pv  has  to  transmit.  In  other  words,  those  vertices  between 
pv  and  v  cannot  directly  reach  v.  Equivalently,  any  vertex  to  the 
right  of  v  can  only  hear  a  message  from  s  through  a  relay  by  pv 
or  a  vertex  to  the  left  of  pv . 

The  NBI  algorithm  is  then  carried  out  in  two  steps.  In  the 
first  step,  the  predecessors  of  certain  vertices  are  obtained  one 
by  one  starting  from  t  moving  toward  s.  Specifically,  the  prede¬ 
cessor  of  t ,  denoted  by  u\  =  pu  is  first  obtained.  If  xUl  <  xs , 

5 We  can  consider  a  reduction  from  MDS  in  graphs  with  a  maximum  degree 
up  to  9  (see  Appendix  D). 


Fig.  3.  1-D  network  (circles  represent  maximum  ranges  under  a  disk  propaga¬ 
tion  model). 


then  the  first  step  terminates.  Otherwise,  the  predecessor  of  u\, 
denoted  by  u2  =  pUl,  is  obtained  and  its  location  compared 
to  xs.  The  same  procedure  continues  until  the  currently  ob¬ 
tained  predecessor  is  to  the  left  of  s  or  is  s  itself.  The  first 
step  thus  produces  a  sequence  of  vertices  m,  u2, . . . ,  ui  with 
ui  =  pt,u2  =  pUl,.-.,ui  =  pUl_1  and  xUl  <  xs.  Then, 
L\  :=  {ui,  ui-i, . . . ,  ui,  t}  is  a  valid  path  from  ui  to  t.  If 
m  =  s,  the  algorithm  terminates,  and  the  thinnest  path  from  s  to 
t  is  given  by  L\ .  Otherwise,  we  carry  out  Step  2  of  the  algorithm 
where  we  find  a  path  from  s  to  ui .  Specifically,  let  V'  denote  the 
set  of  vertices  located  between  ui  and  u/_i  including  ui  but  not 
ui- 1.  Let  E'  denote  the  set  of  all  hyperedges  whose  source  and 
destination  vertices  are  in  V' .  As  shown  in  Appendix  E  on  the 
correctness  of  the  algorithm,  any  hyperpath  L2  from  s  to  ui  in 
the  subhypergraph  H'  =  (V',Er)  concatenated  with  L\  gives  a 
thinnest  path  from  s  tot.  Finding  such  an  L2  can  be  easily  done 
by  a  breadth-first  search  (BFS)  in  H1 .  However,  the  resulting 
time  complexity  is  0(n2).  Hence,  we  propose  a  special  BFS 
procedure  that  reduces  the  time  complexity  to  0(n).  The  trick 
here  is  to  set  up  two  pointers,  ki  and  kr ,  to  the  locations  of  the 
leftmost  and  the  rightmost  vertices  in  V'  that  have  been  discov¬ 
ered.  Due  to  the  geometric  structure  of  the  1  -D  network,  each 
time  we  only  need  to  search  vertices  to  the  left  of  ki  and  vertices 
to  the  right  of  kr.  The  detailed  algorithm  is  given  below. 

1.  Enqueue  s,  set  ki  and  kr  to  the  index  of  s. 

2.  Repeat  until  the  queue  is  empty  or  ui  is  found: 

—  Dequeue  a  vertex  v  and  examine  it 

—  If  v  =  ui,  go  to  step  4. 

—  Otherwise, 

*  While  v  can  reach  Vkr+i 

•  Enqueue  Vkr+i  and  kr  <—  kr  +  1 

•  Set  the  parent  of  Vkr+i  to  v 

*  While  v  can  reach  Vkt  - 1 

•  Enqueue  vj^-i  and  ki  <—  ki  —  1 

•  Set  the  parent  of  v^-i  to  v 

3.  If  the  Queue  is  empty,  return  “no  path  from  s  to  t”. 

4.  Trace  back  to  s  and  return  L2. 

The  following  theorem  establishes  the  correctness  of  the  pro¬ 
posed  NBI  algorithm.  Furthermore,  it  reveals  a  strong  property 
of  the  path  obtained  by  NBI  under  a  disk  propagation  model. 
Specifically,  under  a  disk  propagation  model,  we  define  the  cov¬ 
ered  area  A(L )  of  a  hyperpath  L  =  {ei, . . . ,  em}  as 

m 

A{L)  =  (J  -DSe.,re.  (2) 

i=  1 

where  rCi  is  the  minimum  transmission  range  that  induces  hy¬ 
peredge  ei,  i.e., 

rei  =  max{d(sei,v)}.  (3) 

uETe. 
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Theorem  5  shows  that  the  covered  area  of  the  path  obtained  by 
NBI  is  a  subset  of  the  covered  area  of  any  feasible  path  from  s 
to  t. 

Theorem  5:  NBI  algorithm  finds  the  thinnest  path  L*. 
Furthermore,  under  a  disk  propagation  model,  given  any  valid 
path  L  from  s  to  t ,  we  have  Al*  C  Al . 

Proof:  See  Appendix  E.  ■ 

Theorem  6:  The  time  complexity  of  the  NBI  algorithm  is 
0(n). 

Proof:  The  0(n)  complexity  of  the  first  step  of  NBI  is 
readily  seen.  In  the  second  step,  the  time  complexity  is  domi¬ 
nated  by  updating  the  queue  at  each  iteration.  Let  k  denote  the 
number  of  iterations  in  step  2.  Note  that  we  only  check  rrii  +  2 
vertices  at  iteration  z,  where  rrii  is  the  number  of  new  vertices 
that  have  been  enqueued  at  this  iteration  and  Yli=i  mi  <  in 
Also  k  is  bounded  by  \V'\.  Hence,  the  total  time  complexity  of 
this  step  is  bounded  by  Yli=i  (m*  +  2)  <  3\V'\.  We  thus  arrive 
at  the  theorem.  ■ 


Fig.  4.  Example  where  SPBA  outperforms  TSBA.  There  are  two  paths  from 
s  to  t.  One  goes  through  all  solid  black  hyperedges  to  v  and  then  to  t,  and  the 
other  contains  all  dashed  hyperedges.  The  first  one  is  the  thinnest  path  since  it 
only  covers  six  vertices,  while  the  second  one  covers  all  eight  vertices.  SPBA 
returns  the  first  path  since  its  length  is  8,  while  the  second  one  is  10.  However, 
TSBA  returns  the  second  path  because  the  path  from  s  to  v  is  chosen  to  be  the 
dashed  hyperedge  one  and  is  used  to  generate  the  path  from  s  to  t. 


2(1  + 
with 


2 a)d -approximation  for  d-dimensional  ring  hypergraphs 

 max,.Gy  Ri 


max{minViEv  rr,  minUiX,ev  d(u,v)}  ’ 


Additionally,  the 


ratio  Jf:  of  the  SPBA  algorithms  is  asymptotically  tight  even 


in  2-D  disk  hypergraphs. 
Proof:  See  Appendix  F. 


B.  1 . 5-D  Networks 

We  now  consider  the  1 .5-D  problem  where  in-network  nodes 
are  located  on  a  line  and  eavesdroppers  are  located  in  a  d-di¬ 
mensional  space  that  contains  the  line  network.  We  focus  on  the 
disk  propagation  model.  We  assume  a  unit  cost  for  each  in-net¬ 
work  node  that  hears  the  message  and  a  nonnegative  cost  c  for 
each  eavesdropper  that  hears  the  message,  where  c  can  take  any 
nonnegative  value,  thus  allowing  us  to  model  more  general  sce¬ 
narios  where  overhearing  by  eavesdroppers  can  be  more  costly. 
The  objective  is  to  find  a  path  L*  from  s  to  t  with  the  minimum 
total  cost 

L*  =  arg  min  <  civ)  >  (4) 

L={ei,...,em}  j 

where  c(v)  is  the  cost  for  vertex  v,  and  A(L)  is  the  covered  area 
of  path  L  as  defined  in  (2). 

Based  on  Theorem  5,  the  path  provided  by  NBI  covers  only 
those  essential  vertices  that  must  be  covered  by  any  valid  path.  It 
thus  follows  that  NBI  provides  the  optimal  solution  to  the  1 .5-D 
thinnest  path  problem  without  knowledge  of  the  eavesdroppers 
locations.  More  specifically,  no  algorithm,  even  with  complete 
knowledge  of  the  locations  of  the  eavesdroppers,  can  obtain  a 
thinner  path  than  NBI,  which  does  not  require  location  knowl¬ 
edge  of  the  eavesdroppers. 

V.  Approximation  Algorithms 

In  this  section,  we  introduce  two  approximation  algorithms 
for  the  thinnest  path  problem  and  analyze  their  performance  in 
different  types  of  hypergraphs. 

A.  Shortest-Path-Based  Approximation  Algorithm 

Given  a  general  directed  hypergraph  H  with  source  vertex  s 
and  destination  vertex  t ,  we  set  the  weight  of  a  hyperedge  to  be 
the  number  of  destination  vertices  in  this  hyperedge 

w{e)  =  \Te\.  (5) 

The  shortest  hyperpath  algorithm  from  s  to  t  is  then  obtained 
under  this  weight  definition  as  an  approximation  of  the  thinnest 
path.  The  following  theorem  quantifies  the  performance  of  this 
shortest-path-based  algorithm  (SPBA). 

Theorem  7:  The  SPBA  algorithm  provides  a  y^-ap- 
proximation  for  TP  in  general  directed  hypergraphs,  a 


B.  Tree- Structure-Based  Approximation  Algorithm 

Approximation  occurs  in  two  places  in  SPBA.  First,  the  width 
of  a  path  is  approximated  by  the  sum  of  the  widths  of  the  hy¬ 
peredges  on  that  path.  Second,  the  thinnest  path  to  a  vertex  is 
assumed  to  go  through  the  thinnest  path  to  one  of  its  incoming 
neighbors.  The  first  approximation  can  be  avoided  while  main¬ 
taining  the  polynomial  nature  of  the  approximation  algorithm. 
In  particular,  we  can  ensure  that  the  width  of  a  path  is  correctly 
obtained  by  using  the  set  union  operation  instead  of  summation. 
The  assumption  on  the  tree  structure  of  the  thinnest  paths  allows 
us  to  use  Dijkstra's  algorithm  with  some  modifications.  Specifi¬ 
cally,  for  each  vertex,  we  need  to  store  the  current  thinnest  path 
from  s  to  this  vertex  rather  than  only  the  width  of  this  path 
and  the  parent  of  this  vertex  on  this  path.  This  allows  us  to 
take  the  set  union  operation  when  we  update  the  neighbors  of 
this  vertex.  Given  below  is  the  performance  of  this  tree-struc- 
ture-based  algorithm  (TSBA). 

Theorem  8:  The  TSBA  algorithm  provides  a  ^=== -ap¬ 
proximation  for  general  directed  hypergraphs,  2(1  +  2ai)d-ap- 
proximation  for  d-dimensional  ring  hypergraphs  with 
a  =  j—. max^EV  R* — — TT.  Additionally,  the  ratio 

max{mm7;i  e  y  fi ,  mmU)7;e  y  d{u,v)f  *7’ 

of  the  TSBA  algorithm  is  tight  in  general  directed 
hypergraphs  and  asymptotically  tight  in  disk  hypergraphs  in 
the  worst  case. 

Proof:  See  Appendix  G.  ■ 


C.  Performance  Comparison 

Since  both  SPBA  and  TSBA  are  based  on  a  dynamic  program 
similar  to  the  Dijkstra’s  algorithm  for  shortest  path  in  graphs,  it 
is  not  difficult  to  show  that  the  time  complexities  of  both  algo¬ 
rithms  are  0(<F  +  rzlogrz)  where  =  YheeE  l^el-  Thus,  their 
complexities  are  linear  with  the  size  of  the  given  hypergraph 
H  =  (V,  E),  which  is  order-optimal. 

While  the  approximation  ratio  of  TSBA  is  better  than  that  of 
SPBA,  these  are  worst-case  performances  and  do  not  imply  that 
TSBA  outperforms  SPBA  in  every  case  as  shown  in  Fig.  4. 

Fig.  5  shows  the  average  performance  of  these  two  algo¬ 
rithms.  We  see  that  both  algorithms  have  relatively  small  ap¬ 
proximation  ratios  growing  sublinearly  with  the  number  of  ver¬ 
tices.  In  general,  TSBA  outperforms  SPBA  on  average,  as  also 
demonstrated  in  a  number  of  other  simulation  results  (omitted 
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Number  of  vertices 

Fig.  5.  Average  performance  of  SPBA  and  TSBA  (a  2-D  network  with  n  ver¬ 
tices  uniformly  and  randomly  distributed  on  an  —  x  —  square  with  p  =  1.5;  the 
maximum  range  of  each  vertex  is  randomly  chosen  from  interval  [Rm  in ,  i?max] 
with  Rmin  =  1,  f?max  =  5;  average  taken  over  1000  such  random  2-D  disk 
hypergraphs). 

due  to  the  space  limit).  However,  the  performance  of  SPBA  has 
a  smaller  variance  than  that  of  TSBA.  This  is  mainly  due  to  the 
fact  that  the  thinnest  path  itself  has  a  larger  variance  than  the 
shortest  path  (as  confirmed  in  our  simulations),  and  TSBA  often 
returns  the  thinnest  path  rather  than  the  shortest  path. 

VI.  Conclusion 

We  studied  the  complexity  and  developed  optimal  and  ap¬ 
proximation  algorithms  for  the  thinnest  path  problem  for  secure 
communications  in  wireless  ad  hoc  networks.  In  establishing 
the  NP-completeness  of  the  problem,  our  techniques  of  using 
directed  crosses  and  exposed  disk  hypergraphs  may  spark  new 
tools  for  complexity  studies  in  geometrical  hypergraphs  and 
graphs.  The  bounding  techniques  and  the  use  of  sphere  packing 
results  in  analyzing  the  performance  of  the  two  approximation 
algorithms  may  also  find  other  applications  in  algorithmic  anal¬ 
ysis.  Whether  the  proposed  approximation  algorithm  TSBA  of¬ 
fers  the  optimal  approximation  ratio  is  still  an  open  question  that 
requires  further  investigation. 

Appendix  A 
Proof  of  Theorem  1 

A.  Reduction  From  MDS  to  TP  in  a  General  Directed 
Hypergraph  H\ 

Consider  the  MDS  problem  in  a  graph  G  with  n  vertices 
vi , . . . ,  vn.  We  construct  a  directed  hypergraph  Hi  based  on  G 
as  follows.  The  vertex  set  of  Hi  includes  the  n  vertices  of  G 
augmented  by  a  destination  vertex  un+i  and  n  super  vertices 
v{ , . . . ,  v!^ .  A  super  vertex  vf  corresponds  to  the  normal  vertex 
Vi  and  is  a  set  of  ns  normal  vertices.  The  hyperedges  in  Hi  are 
all  rooted  at  the  normal  vertices  vi,...,vn.  Specifically,  rooted 
at  Vi  (1  <  i  <  n)  are  ki  + 1  directed  hyperedges,  where  ki  is  the 
degree  of  Vi  in  G.  Each  hyperedge  rooted  at  Vi  has  two  destina¬ 
tions:  Vi+i  and  a  super  vertex  v*j  whose  corresponding  normal 
vertex  Vj  dominates6  Vi  in  the  original  graph  G.  Fig.  6  is  an  ex¬ 
ample  illustrating  the  construction  of  Hi  from  G. 

From  the  construction  of  Hi ,  we  see  that  any  path  from  vi 
to  un+i  must  traverse  through  all  normal  vertices  one  by  one. 
There  are  multiple  hyperedges  leading  from  Vi  to  V{+i,  each 
involving  a  super  vertex  that  corresponds  to  a  dominating  node 
of  Vi  in  G.  Thus,  choosing  a  hyperedge  going  from  Vi  to  vi+\ 
is  equivalent  to  choosing  a  dominating  node  of  Vi  in  G.  Since 
every  path  from  vi  to  un+i  includes  all  the  n+ 1  normal  vertices, 

6  A  vertex  in  a  graph  is  dominated  by  itself  and  any  of  its  one-hop  neighbors. 


Fig.  6.  Construction  of  ifi  from  G :  (a)  the  graph  G;  (b)  the  hypergraph  Hi 
(vi  is  dominated  by  tq  and  v3  in  G.  We  thus  have  two  hyperedges  rooted  at  tq 
in  Hi :  One  reaches  ( v2 ,  v* ),  the  other  (u2 ,  nf ) .). 

the  thinnest  path  is  given  by  the  one  with  the  minimum  number 
of  super  vertices,  thus  leading  to  the  MDS  in  G.  At  this  point, 
the  size  ns  of  a  super  vertex  can  be  any  positive  integer.  As  will 
become  clear  later,  to  implement  Hi  under  a  2-D  disk  model, 
additional  normal  vertices  need  to  be  added.  As  a  consequence, 
paths  from  vi  to  un+i  may  include  different  numbers  of  normal 
vertices.  To  preserve  the  reduction,  we  need  to  make  sure  that 
the  width  of  a  path  is  dominated  by  the  number  of  super  vertices 
it  covers.  This  can  be  achieved  by  choosing  an  ns  sufficiently 
large  (see  Appendix  A-D). 

The  following  lemma  formally  establishes  the  correctness  of 
the  reduction. 

Lemma  3:  There  is  a  dominating  set  with  size  k  in  G  if  and 
only  if  there  is  a  path  from  vi  to  un+i  in  Hi  with  width  kns  + 
n  +  1. 

Proof:  First,  assume  that  G  has  a  dominating  set  S  with 
size  k.  By  the  definition  of  dominating  set,  for  each  vertex  Vi 
in  G ,  there  is  a  vertex  Vj  e  S  that  dominates  Vi.  From  the  con¬ 
struction  of  Hi ,  there  exists  a  hyperedge  e*  (i  =  1 , . . . ,  n)  in 
Hi  directed  from  Vi  to  vertex  V{+i  and  super  vertex  vs-  corre¬ 
sponding  to  the  dominating  node  v3  in  G.  Thus,  the  hyperpath 
{ei, . . . ,  en}  is  a  path  from  vi  to  un+i  with  width  kns  +  n  +  1. 
The  width  comes  from  the  fact  that  all  n  +  1  vertices  in  Vr1  are 
on  the  path  along  with  k  super  vertices,  each  consisting  of  ns 
normal  vertices. 

Conversely,  assume  that  there  exists  a  path  from  vi  to  un+i 
in  Hi  with  width  kns  +  n  +  1.  Based  on  the  construction  of 
H ,  every  path  from  vi  to  un+i  consists  of  n  hyperedges  rooted 
at  each  of  the  n  normal  vertices  vi, ...  ,vn.  Thus,  a  path  with 
width  kns  +  n  +  1  must  contain  k  super  vertices.  From  the 
construction  of  the  hyperedges,  we  conclude  that  the  vertices  in 
G  that  correspond  to  those  k  super  vertices  along  the  given  path 
form  a  dominating  set  with  size  k.  ■ 

B.  2-D  Grid  Representation  of  Hi 

The  directed  hypergraph  Hi  obtained  above  does  not  sat¬ 
isfy  the  geometric  properties  of  2-D  disk  hypergraphs  (see 
Section  II).  To  prove  Theorem  1,  we  need  to  modify  Hi  to  a 
2-D  disk  hypergraph  H2  while  preserving  the  reduction  from 
MDS  in  G.  Our  approach  is  to  realize  the  topological  structure 
of  each  hyperedge  in  Hi  by  adding  additional  vertices  with 
carefully  chosen  locations  and  maximum  ranges  to  lead  from 
the  source  vertex  to  the  destination  vertices  of  this  hyperedge. 
The  number  of  additional  vertices,  however,  should  be  kept  at 
a  polynomial  order  with  the  problem  size  to  ensure  the  polyno¬ 
mial  nature  of  the  reduction.  This  can  be  achieved  by  adding 
vertices  on  a  2-D  grid  with  a  constant  grid  spacing,  which 
allows  a  constant  maximum  range,  thus  polynomially  many 
additional  vertices.  The  detailed  implementation  of  Hi  under 
a  2-D  disk  model  is  given  in  Appendix  A-C.  As  a  preparatory 
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Fig.  7.  2-D  grid  representation  of  (the  two  hyperedges  rooted  at  v\  from 
the  example  given  in  Fig.  6  are  illustrated  in  green  and  blue,  respectively). 


step,  we  show  in  this  section  that  the  hyperedges  in  Hi  can  be 
represented  by  line  segments  of  a  2-D  grid  with  a  constant  grid 
spacing. 

We  first  embed  the  normal  vertices  of  Hi  evenly  in  a  hori¬ 
zontal  line  in  a  2-D  space  (see  Fig.  7  for  an  illustration).  Below 
this  line  is  a  2 n2  x  4 n2  unit  grid.  There  are  4 n  vertical  lines 
between  V{  and  i(l  <  i  <  n)  that  are  partitioned  into  three 
zones  (Cf,  Cf,  Cf )  of  n,  2 n,  and  n  vertical  lines,  respectively. 
The  super  vertices  are  embedded  evenly  on  a  horizontal  line 
below  the  grid.  The  horizontal  position  of  super  vertex  vf  is  be¬ 
tween  Vi  and  Vi+i. 

Next,  we  specify  how  a  hyperedge  traverses  the  grid  from  its 
source  vertex  to  its  destination  vertices.  Recall  that  every  hy¬ 
peredge  in  Hi  is  directed  from  a  normal  vertex  v%  to  a  super 
vertex  vf  and  the  next  normal  vertex  vi+i.  To  preserve  the  re¬ 
duction,  we  need  to  ensure  that  each  hyperedge  can  only  reach 
its  normal  vertex  destination  after  reaching  its  super  vertex  des¬ 
tination.  To  facilitate  the  implementation  around  the  super  ver¬ 
tices  (see  Appendix  A-C.2),  we  designate  the  middle  zone  Cf 
between  Vi  and  for  traveling  down  to  super  vertex  vf  and 
then  up  to  the  corresponding  normal  vertex  destination  (see  re¬ 
gion  Cf  in  Fig.  7).  Each  hyperedge  involving  vf  has  two  desig¬ 
nated  vertical  lines  in  Cf  (one  for  going  down  to,  the  other  going 
up  from,  the  super  vertex).  To  connect  the  designated  vertical 
lines  in  zone  Cf  with  the  source  vertex  and  then  to  the  normal 
destination  vertex,  we  designate  two  horizontal  lines  for  each 
hyperedge.  The  traverse  of  the  hyperedge  is  completed  by  des¬ 
ignating  one  vertical  line  in  C\  and  one  in  Cf  to  connect  the 
normal  vertices  with  the  corresponding  designated  horizontal 
lines.  Since  there  are  at  most  n2  hyperedges,  the  designed  grid 
size  is  sufficient  to  ensure  that  each  hyperedge  traverses  through 
a  distinct  set  of  line  segments  in  the  grid. 

C.  Implementing  Hi  Under  a  2-D  Disk  Model 

Based  on  the  2-D  grid  representation  of  Hi,  we  can  construct 
a  2-D  disk  hypergraph  H2  that  preserves  the  reduction.  Specif¬ 
ically,  we  place  a  sequence  of  evenly  spaced  normal  vertices 
with  a  constant  maximum  range  along  the  line  segments  in  the 
grid  that  form  each  hyperedge  of  Hi.  The  distance  between  two 
adjacent  vertices  is  set  to  their  maximum  range.  The  constant 
maximum  range  can  be  set  sufficiently  small  (say,  |)  to  avoid 
overhearing  across  vertices  on  different  hyperedges  that  may 
render  the  reduction  invalid.  There  are  two  issues  that  remain 
to  be  addressed:  the  implementation  of  crosses  and  that  around 
super  vertices. 


Fig.  8.  Disk  hypergraph  implementation  of  a  directed  cross  where  the  circles 
represent  the  maximum  range  of  vertices  (messages  transmitted  on  the  blue  line 
can  be  heard  by  nodes  on  the  red  line,  but  not  vice  versa). 

1)  Implementation  of  Crosses:  The  line  crossing  in  the  grid 
representation  of  Hi  makes  overhearing  across  hyperedges  in¬ 
evitable.  However,  by  exploiting  the  freedom  of  choosing  the 
maximum  range  for  each  vertex,  we  can  implement  directed 
crosses  that  allow  us  to  preserve  the  reduction.  Specifically, 
when  two  line  segments  in  the  grid  representation  cross,  we  can 
choose  the  maximum  ranges  of  the  vertices  along  these  two  lines 
in  such  a  way  that  messages  transmitted  over  one  line  can  be 
heard  by  vertices  on  the  other  but  not  vice  versa.  A  specific  im¬ 
plementation  is  given  in  Fig.  8. 

Next,  we  show  how  carefully  choosing  the  direction  of  each 
cross  allows  us  to  preserve  the  reduction.  The  cross  directions 
are  defined  by  assigning  a  level  index  to  each  line  segment  in  the 
grid  representation.  Specifically,  for  a  hyperedge  rooted  at  Vi  in 
Hi ,  its  line  segments  before  and  after  reaching  the  super  vertex 
destination  have  levels  i  andz+1,  respectively.  Then,  each  cross 
has  a  direction  pointing  from  the  higher  level  segment  to  the 
lower  one  (i.e.,  messages  transmitted  on  the  higher-level  seg¬ 
ment  can  be  heard  by  the  vertices  along  the  lower-level  seg¬ 
ment,  but  not  vice  versa).  If  the  two  segments  have  the  same 
level,  the  direction  of  the  cross  can  be  arbitrary.  To  see  that  this 
directed  implementation  of  crosses  preserves  the  reduction,  we 
only  need  to  notice  that  any  path  from  vi  to  un+i  still  needs  to 
go  through  all  the  n  normal  vertices  one  by  one  and  must  reach 
a  super  vertex  before  reaching  the  next  normal  vertex. 

2)  Implementation  Around  Super  Vertices:  Recall  that  a 
super  vertex  in  Hi  is  a  set  of  ns  normal  vertices  that  have  no 
outgoing  hyperedges.  It  can  be  implemented  by  ns  points  with 
zero  maximum  range  and  located  sufficiently  close  to  each 
other  (so  that  any  path  from  vi  to  un+i  in  H2  includes  either 
all  of  them  or  none  of  them). 

Consider  first  the  implementation  of  one  incoming  hyperedge 
to  a  super  vertex  vf  Recall  that  in  the  2-D  grid  representa¬ 
tion  of  Hi ,  a  hyperedge  approaches  and  leaves  vs-  through  two 
vertical  lines  in  zone  Cf  (see  Fig.  7).  One  implementation  of 
this  U-turn  around  vf  is  to  add  six  normal  vertices  with  spe¬ 
cific  maximum  ranges  and  locations.  As  shown  in  Fig.  9,  these 
six  vertices  include  three  anchor  vertices  vf ,  and  u f  with 
maximum  range  r,  two  interface  vertices  vf  and  that  con¬ 
nect  with  the  grid,  and  a  bridging  vertex  /in,  all  with  maximum 
range  | .  The  value  of  r  and  the  connection  with  the  grid  will  be 
specified  later. 

A  challenge  remains  in  the  implementation  of  up  to  n  in¬ 
coming  hyperedges  to  the  same  super  vertex.  Note  that  under 
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Fig.  9.  Implementation  of  one  hyperedge  passing  through  a  super  vertex. 
Starting  from  ,  the  message  traverses  to  t/f  through  /in ,  uf ,  u\ ,  .  The 

super  vertex  hears  the  message  in  the  transmission  from  u\  to  u  f . 


Fig.  10.  Implementation  of  the  second  incoming  hyperedge  to  a  super  vertex. 

a  2-D  disk  model,  one  can  at  most  have  five  vertices  (even  with 
arbitrary  ranges)  that  reach  a  common  sixth  vertex  but  not  each 
other.  The  key  to  circumvent  this  difficulty  is  to  allow  directed 
overhearing,  similar  to  the  idea  behind  the  implementation  of 
the  crosses.  Specifically,  the  reduction  is  preserved  as  long  as  a 
hyperedge  rooted  at  Vj  cannot  overhear  a  message  transmitted 
over  a  hyperedge  rooted  at  V{  for  any  i  <  j.  The  detailed  imple¬ 
mentation  is  as  follows.  The  fist  step  is  to  designate  the  vertical 
lines  in  zone  C?  to  the  incoming  hyperedges  of  vs-  based  on 
the  indices  of  their  source  vertices.  Specifically,  the  incoming 
hyperedge  with  the  smallest  source  vertex  index  takes  the  two 
centermost  lines  in  Cj,  and  so  on.  Consider  first  the  implemen¬ 
tation  of  the  two  incoming  hyperedges  (say,  e\  and  e2)  with  the 
smallest  source  vertex  indices.  As  shown  in  Fig.  10,  we  first 
implement  e\  as  described  above  (see  Fig.  9).  The  structure  of 
the  implementation  of  e2  is  similar  except  that  the  maximum 
range  of  the  anchor  vertices  u2  ,  u2,  and  u  2  is  set  to  4 r  to  pre¬ 
vent  unwanted  overhearing.  As  a  consequence,  more  bridging 
vertices  (^21,^22,^23  with  maximum  range  r,  and  2 r,  re¬ 
spectively)  are  needed  to  connect  the  interface  vertex  v2  to  the 
anchor  vertex  u2  .  Note  that  no  vertices  along  e2  (the  centers  of 
the  blue  circles  in  Fig.  10)  are  in  the  range  of  any  vertices  along 
ei  (the  green  circles).  The  correct  direction  of  overhearing  is 
thus  ensured. 

The  same  procedure  continues  for  any  additional  incoming 
hyperedges  to  vf ,  in  the  ascending  order  of  their  source  vertex 
indices  in  Hi.  Note  that  the  range  of  the  anchor  vertices  in  the 
kth  hyperedge  is  4fcr,  growing  exponentially  with  k.  The  max¬ 
imum  ranges  (specifically,  r, . . . ,  ^-r)  of  the  bridging  ver¬ 
tices  {fiki}  are  chosen  to  preserve  the  polynomial  nature  of  the 
reduction.  In  this  way,  the  number  of  additional  vertices  for  im¬ 
plementing  the  kth  hyperedge  is  2k  +  4,  and  the  total  number  of 
additional  vertices  around  one  super  vertex  is  at  most  n 2  +  3 n. 

Next,  we  consider  the  value  of  r,  which  should  be  set  suf¬ 
ficiently  small  to  avoid  overhearing  across  hyperedges  leading 
to  different  super  vertices.  Note  that  the  width  of  the  area  cov¬ 
ered  by  the  additional  vertices  around  a  super  vertex  is  4  times 
the  largest  maximum  range  of  the  anchor  vertices.  We  thus  set 
r  =  4-nn,  considering  the  distance  between  two  adjacent  super 
vertices  being  An. 
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Fig.  11.  Consider  first  the  downward  part  from  the  grid  to  a  left  interface  vertex 
i/T .  Let  O  i  denote  the  location  of  the  last  vertex  on  the  designated  vertical 
line  in  the  grid,  and  02  the  location  of  i/E .  The  circles  centered  at  0\  and  02 
represent  their  maximum  ranges.  Let  A,  B  and  C,  D  denote  the  intersecting 
points  of  these  two  circles  with  the  horizontal  lines  at  their  centers.  Let  Ei 
denote  the  intersection  between  circle  O i  and  line  Cfi  02 .  Let  di  and  d[  denote 
the  distance  between  E\  and  the  two  lines  AC  and  BD,  respectively.  Next,  we 
draw  a  circle  with  radius  r±  =  min{di ,  d[  }  centered  at  Ex .  Let  E2  denote  the 
intersection  between  circle  E\  and  line  0i02,  and  a  similar  circle  centered  at 
E2  is  drawn.  This  procedure  is  repeated  to  generate  a  sequence  of  circles  until 
the  last  generated  circle  covers  02.  This  sequence  of  circles  { E j}  gives  the 
locations  and  the  maximum  ranges  of  the  vertices  connecting  the  grid  and  . 
The  upward  part  from  vf  to  the  grid  is  done  with  the  same  procedure  except 
starting  from  vf . 

The  last  issue  is  to  connect  the  interface  vertices  with  the  grid. 
Each  interface  vertex  needs  to  be  connected  with  a  designated 
vertical  line  in  Cj.  While  the  vertical  lines  in  Cj  are  evenly 
spaced,  the  horizontal  positions  of  the  interface  vertices  have  an 
exponential  structure  due  to  the  exponentially  growing  range  of 
the  anchor  vertices.  Furthermore,  the  vertices  realizing  the  ver¬ 
tical  lines  in  the  grid  have  a  constant  range,  whereas  the  inter¬ 
face  vertices  have  an  exponentially  smaller  range  of  r  =  A~nn. 
If  we  connect  them  using  a  sequence  of  vertices  with  a  constant 
range,  unwanted  overhearing  will  occur  near  the  interface  ver¬ 
tices.  On  the  other  hand,  connecting  them  using  vertices  with 
range  r  results  in  an  exponential  number  of  additional  vertices. 
To  preserve  the  correctness  and  the  polynomial  nature  of  the  re¬ 
duction,  we  propose  the  scheme  detailed  in  Fig.  11. 

Since  the  generated  sequence  of  circles  {Ei}  are  within  the 
boundary  given  by  lines  AC  and  BD  and  the  boundary  lines 
corresponding  to  different  interface  vertices  do  not  cross  (see 
Fig.  12),  the  above  scheme  does  not  introduce  overhearing,  thus 
preserving  the  reduction.  The  polynomial  nature  of  the  reduc¬ 
tion  can  be  shown  based  on  the  following  lemma. 

Lemma  4:  Consider  the  geometrical  scheme  described  in 
Fig.  11.  Assume  LCAB  >  J.  The  number  of  circles  {Ei},  de¬ 
noted  by  k,  satisfies  k  <  2(1°g^1j^°gfi2)  e  +  1  when  Ri  /  R2, 
and  k  <  +  1  when  Ri  =  R2,  where  R\ ,  R2  denote  the 

radii  of  circles  0\  and  02 ,  and  L  the  distance  between  lines 
AB  and  CD. 

Proof:  Assume  first  Ri  f  R2 .  Without  loss  of  generality, 
assume  Ri  >  R2.  Since  lines  AB  and  CD  are  parallel,  the  three 
lines  AC,  0\02 ,  and  BD  intersect  at  one  point,  denoted  by  O 
in  Fig.  11.  Let  a ,  (3,  and  9  denote  the  angles  ZOiOA,  LOAB, 
and  LOO\B,  respectively. 
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Fig.  12.  Connecting  the  interface  vertices  with  the  grid. 


It  can  be  shown  that  all  the  circles  {Ef\  are  tangential  to  the 
same  boundary  line.  Without  loss  of  generality,  assume  that  the 
tangential  line  is  AC ,  i.e.,  d{  <  d[  and  ri  =  cfy  Based  on  simple 
geometry,  the  lengths  of  the  line  segments  of  {OEi}  form  an 
equal  ratio  sequence 

OEi+i  =  OEi  —  ri  =  OE{(  1  —  sin  ct) 

with  OEi  =  00\  —  R\.  We  thus  have 

OEi+ 1  =  (OOi  -  i?i)(l  -  sin  a)*. 

Based  on  the  stopping  condition  of  the  procedure,  the  number  k 
of  circles  is  given  by  the  minimum  index  i  such  that  OE{+ 1  < 
002-  We  thus  have 


k  =  min{i  G  N  :  0Ei+ 1  <  00 2} 

=  min{i  G  N  :  ( 00\  —  ffy)(l  —  sin  a)*  <  00 2} 

•  i  '  r  M  ■  /  OOi-fli  1 
=  mm  <  1  G  N  :  1  <  - — — - ^ > 

[  log(l  —  sma)  I 

<  \og(002/ 00\) /  log(l  —  sin  a)  +  1. 


(6) 


Since  AOO2D  and  AOOiB  are  similar  triangles,  the  ratio 
equals  the  ratio  Also  because  —  log(l  —  x)  >  x  for 
0  <  x  <  1,  (6)  can  be  written  as 


k  <  (logfth  —  logi?2)/  sinct  +  1.  (7) 

Because  0\  AO  is  a  triangle  and  (3  >  J,  the  value  of  sin  a  can 
be  lower-bounded  as  follows: 


sin  a  = 


Ri 

OOi 


sin  (3  > 


Ri  V2 

OOi  2  - 


Ri  —  i?2  v^2 

0i02  2“' 


(8) 


Furthermore,  since  0  =  a  +  /?  >  (3,  the  length  of  O1O2  = 
has  an  upper  bound:  O1O2  <  <  \/2T.  Hence,  (8) 

leads  to 


sin  a  >  (R1-R2)/{2L).  (9) 

Substituting  (9)  into  (7),  we  have 


k  <  2L(logi^i  -  log  tf2)/(#i  -  ^2)  +  1. 

Consider  next  =  i?2-  The  sequence  of  circles  {^}  have 
the  same  radius  Ri.  Since  O1O2  =  <  2L,  the  bound 

k  <  ^  +  1  holds.  ■ 

To  satisfy  the  assumption  of  LCAB  >  J  in  Lemma  4,  we 
set  the  distance  between  the  last  horizontal  line  of  the  grid  and 
the  horizontal  line  of  super  vertices  to  n.  This  ensures  that  angle 


LBAO  <  Note  that  in  the  downward  part  from  the  grid  to 
a  left  interface  vertex  i/f,  ffy  is  a  constant  and  R2  =  4  . 

Hence,  the  bound  on  k  given  in  Lemma  4  can  be  written  as 

2(log  Ri  +  n  log  4  —  log  n) 

k  <  - — - 7 - — - L  +  1 

—  4  nn/2 

2  (log  +  nlog4) 

-  n  +  1 

which  is  in  the  order  of  0{n2) .  A  similar  argument  can  be  made 
for  the  upward  part  where  R\  =4 ~nn  and  R2  is  a  constant.  The 
same  holds  for  ffy  =  #2 .  Hence,  the  total  number  of  additional 
vertices  to  connect  the  grid  to  the  interface  vertices  of  a  super 
vertex  is  in  the  order  of  0(n3). 

D.  Reduction  From  MDS  to  TP  in  the  2-D  Disk  Hypergraph 

H2 

With  H2  constructed,  we  now  establish  the  correctness  of  the 
reduction  from  the  MDS  in  G  to  the  TP  from  v\  to  un+i  in  H2. 

Lemma  5:  Let  ns  =  n 2  +  1  where  n 2  is  the  total  number  of 
normal  vertices  in  H2 .  There  is  a  dominating  set  with  size  k  in 
G  if  and  only  if  there  is  a  path  from  v\  to  un+i  in  H2  with  width 
between  kns  and  ( k  +  l)ns  —  1. 

Proof:  The  chosen  value  of  ns  ensures  that  the  width  of  a 
path  from  v\  to  un+i  is  dominated  by  the  number  of  super  ver¬ 
tices  that  it  covers.  The  correctness  of  the  reduction  thus  follows 
from  the  same  arguments  in  the  proof  of  Lemma  3  based  on  the 
construction  of  H2 .  ■ 

The  polynomial  nature  of  the  reduction  is  clear  from  the  con¬ 
struction  of  H2 .  We  thus  arrive  at  Theorem  1 . 

Appendix  B 
Proof  of  Lemma  1 

Consider  a  TP  problem  from  s  to  t  in  a  k- D  exposed  disk 
hypergraphs  H  =  (V,  E ).  We  construct  a  k- D  UDH  H'  as  fol¬ 
lows.  First,  the  normal  vertex  set  V'  of  H'  is  given  by  V,  ex¬ 
cept  that  the  ranges  of  any  v'  G  V'  equals  maxvGy  Rv  Next, 
for  each  vertex  v'  G  V\  we  place  a  super  vertex  in  <3^  (i.e.,  the 
exposed  area  of  the  corresponding  vertex  in  H)  that  contains 
\V\  +  1  normal  vertices  located  sufficiently7  close  to  each  other. 
The  super  vertices  have  the  same  range  as  the  normal  vertices 
in  V',  ensuring  H'  is  a  UDH.  The  reduction  can  thus  be  seen  by 
noticing  that  while  the  enlarged  ranges  introduce  additional  hy¬ 
peredges  in  H',  these  hyperedges  cannot  be  on  a  thinnest  path 
due  to  the  fact  that  they  all  contain  at  least  one  super  vertex. 

Appendix  C 
Proof  of  Lemma  2 

In  this  proof,  we  modify  the  2-D  disk  hypergraph  H2  in  the 
proof  of  Theorem  1  to  a  2-D  exposed  disk  hypergraph  iT3  while 
preserving  the  polynomial  reduction.  Based  on  the  definition, 
a  sufficient  condition  for  a  2-D  disk  hypergraph  to  be  exposed 
is  that  none  of  the  maximum  range  disks  are  completely  inside 
any  other.  The  vertices  in  H2  for  realizing  the  line  segments 
of  the  grid  satisfy  this  condition.  We  only  need  to  modify  the 
implementations  of  the  crosses  and  around  the  super  vertices. 

A.  Implementation  of  Crosses 

In  the  implementation  of  directed  crosses  in  H2  (see  Fig.  8), 
some  vertices  on  the  line  with  a  lower-level  index  may  have  an 

7The  |  V  |  + 1  normal  vertices  are  sufficiently  close  such  that  any  transmission 
from  one  of  these  vertices  to  a  vertex  outside  this  super  vertex  reaches  all  the 
|  V  |  +  1  normal  vertices  in  this  super  vertex. 
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BC  =  Rb  =  Rc  =  R  tan  0  ~  0.5543# 

x  N  BD  =  (1  —  j^)R  cz  0.4283R 


Fig.  1 3 .  To  implement  a  directed  cross  shown  in  (a),  we  first  implement  a  vertex 
for  the  blue  line  with  maximum  range  R  at  location  A  (the  blue  circle)  shown  in 
(b).  Next,  we  draw  a  perpendicular  bisector  between  A  and  E  (the  right  inter¬ 
secting  point  of  the  circle  with  the  line).  On  this  vertical  line,  we  find  two  points 
B  and  C  such  that  £BAE  =  L CAE  =  29°.  At  each  point,  we  put  a  vertex 
for  the  red  line  with  radius  equal  to  the  length  of  BC  [illustrated  by  the  two 
red  circles  in  (b)].  Simple  geometry  calculation  leads  to  BD  <  BC  <  AB  = 
BE.  This  ensures  that  vertices  B  and  C  are  exposed  yet  cannot  overhear  ver¬ 
tices  located  at  A  and  E.  We  complete  the  implementation  by  adding  vertices 
on  the  vertical  line  BC  and  the  horizontal  line  AE  [see  (c)  and  (d)].  Note  that 
to  preserve  the  exposure  of  vertices  B  and  C,  the  maximum  ranges  of  vertices 
from  point  E  to  the  right  side  need  to  be  enlarged  gradually  to  the  constant  max¬ 
imum  range  of  normal  vertices  on  the  grid  (this  only  requires  a  constant  number 
of  additional  vertices). 


empty  exposed  area  (see  the  red  disks  in  Fig.  8  that  are  com¬ 
pletely  covered  by  blue  ones).  To  implement  a  direct  cross  in  a 
2-D  exposed  disk  hypergraph,  the  maximum  ranges  of  vertices 
on  the  line  with  a  lower-level  index  need  to  be  small  enough 
to  preserve  the  direction  of  the  cross,  but  also  large  enough  to 
make  the  vertices  exposed.  We  propose  the  scheme  described  in 
Fig.  13. 

B.  Implementation  Around  Super  Vertices 

In  the  previous  implementation  around  a  super  vertex  v? ,  all 
the  vertices  are  exposed  except  the  anchor  vertices  {u^,uf } 
and  the  bridging  vertices  {gik}-  However,  we  notice  that  these 
vertices  would  all  be  exposed  if  there  were  no  interface  vertices. 
Our  solution  is  thus  to  move  all  the  interface  vertices  away  from 
their  original  positions  by  a  constant  distance  and  add  a  constant 
number  of  vertices  to  connect  each  new  interface  vertex  to  the 
bridging  vertex  or  the  anchor  vertex  on  the  right  side.  A  detailed 
implementation  is  shown  in  Fig.  14. 


Fig.  14.  Interface  vertex  on  the  left  side  is  replaced  by  three  vertices  with  max¬ 
imum  ranges  r,  |  r,  and  J,  respectively.  These  three  vertices  are  located  on  a 
vertical  line  to  the  left  side  of  the  original  location  of  the  interface  vertex  with 
a  distance  of  r.  An  interface  vertex  on  the  right  side  is  replaced  by  two  vertices 
with  maximum  range  ^  located  on  a  vertical  line  to  the  right  side  of  the  orig¬ 
inal  location  of  the  interface  vertex  with  distance  ^ .  Under  this  implementation, 
the  exposed  areas  of  the  anchor  and  bridging  vertices  are  right  above  the  point 
where  they  are  tangential  with  the  horizontal  line  of  the  super  vertices  (as  illus¬ 
trated  by  the  arrows). 


=> 


Fig.  15.  An  implementation  of  a  cross  in  3-D  UDG. 


Fig.  16.  Implementation  around  the  super  vertices  in  UDG. 


In  this  implementation,  there  is  no  overhearing  between  vertices 
on  these  two  line  segments  at  all.  Since  G  has  a  maximum  de¬ 
gree  of  3,  there  are  at  most  four  hyperedges  passing  through 
a  super  vertex.  It  can  be  easily  implemented  without  any  un¬ 
wanted  overhearing  (see  Fig.  16).  To  prevent  the  super  vertices 
from  relaying  messages,  we  place  a  mega  vertex  besides  each 
super  vertex.  This  mega  vertex  is  only  within  the  range  of  this 
super  vertex  and  contains  more  normal  vertices  than  the  total 
number  of  normal  vertices  in  the  reduced  graph  (including  the 
normal  vertices  contained  in  all  the  super  vertices  but  not  those 
in  other  mega  vertices).  In  this  way,  a  path  via  any  super  vertex 
covers  at  least  one  mega  vertex,  thus  cannot  be  the  thinnest  path. 
Fig.  16  illustrates  the  implementation  around  a  super  vertex.8 
The  correctness  of  the  reduction  follows  from  the  same  argu¬ 
ments  as  in  the  proof  of  Lemma  3. 

Appendix  E 
Proof  of  Theorem  5 


Appendix  D 
Proof  of  Theorem  4 

Consider  an  MDS  problem  in  a  graph  G  with  a  maximum  de¬ 
gree  of  3.  We  first  follow  the  first  two  steps  in  the  proof  of  The¬ 
orem  1  to  build  the  grid  representation  of  hypergraph  H\.  Note 
that  due  to  the  unit  range  of  all  vertices,  we  set  the  size  of  the 
grid  to  a  constant  greater  than  1  (say,  5)  to  avoid  unwanted  over¬ 
hearing.  Next,  we  implement  this  representation  in  a  3-D  UDG 
while  preserving  the  reduction.  Any  line  segment  of  hyperedges 
in  Hi  is  replaced  by  a  sequence  of  unit  disks,  one  just  touching 
the  another.  Any  cross  between  two  line  segments  can  be  easily 
implemented  by  using  the  third  dimension,  as  shown  in  Fig.  15. 


We  first  show  that  as  long  as  there  exists  a  path  from  s 
to  t,  there  exists  a  path  from  s  to  ui  that  traverses  only  the 
subhypergraph  H' .  This  can  be  shown  by  noticing  that  ui  must 
hear  the  message  from  s  before  u/_i  and  any  vertex  to  the  right 
of  ui-i.  This  is  due  to  the  mono  tonicity  of  wireless  broadcast 
and  the  definition  of  predecessor.  Consequently,  there  must 
exist  a  path  from  s  to  ui  in  H' .  Since  V'  is  covered  by  the 
hyperedge  leading  from  ui  to  in  Li,  the  concatenation  of 
L\  with  any  path  to  ui  in  H'  covers  the  same  set  of  vertices. 

8We  can  consider  reduction  from  MDS  in  graphs  with  a  maximum  degree  up 
to  9.  In  this  case,  there  are  at  most  10  incoming  hyperedges.  Along  with  the 
mega  vertex,  they  can  be  packed  around  a  super  vertex  without  overhearing. 
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Specifically,  the  cover  of  the  path  returned  by  NBI  is  the  set  of 
vertices  located  between  (and  including)  ui  and  t.  Since  any 
path  from  s  to  t  covers  this  set  of  vertices,  the  correctness  of 
the  algorithm  is  established. 

Next,  we  prove  the  property  of  Al*  under  the  disk  propa¬ 
gation  model.  We  first  state  the  following  lemma  that  follows 
directly  from  triangle  inequality. 

Lemma  6:  Let  D\  and  D2  denote  two  closed  balls  in  Rd  with 
radii  r i  and  r2,  respectively.  Let  a  denote  the  distance  between 
the  centers  of  D\  and  D2 .  If  0  <  a  <  \r±  —  r2\,  then  D2  cDi. 

Based  on  Lemma  6,  for  any  vertex  v  between  ui 
and  ui-i,  we  have  DVtRv  C  Therefore, 

Al *  =  ALl  =  ULi  ,<*(«*,«„_!)  (let  M0  =  t).  Next, 
consider  an  arbitrary  path  L  from  s  to  t.  We  show  that  for  any 
uk(k  =  ,:i  C  Al.  Specifically,  since 

Uk- 1  must  first  hear  the  message  from  Uk  or  a  vertex  to  the  left 
of  Uk,  DUkjd(Ukj  Uk_1)  is  a  subset  of  the  covered  area  of  this  hop 
in  L  based  on  Lemma  6.  This  completes  the  proof. 

Appendix  F 
Proof  of  Theorem  7 

A.  For  General  Directed  Hypergraphs 

Let  Li  denote  the  path  from  s  to  t  provided  by  SPBA  and 
Lopt  =  {ei,  e2, . . . ,  ek}  the  thinnest  path.  If  multiple  thinnest 
paths  exist,  let  Lopt  be  the  one  with  the  minimum  number  of 
hyperedges.  Let  C(L)  denote  the  length  (i.e.,  the  sum  of  hyper¬ 
edge  weights)  of  L. 

Since  each  vertex  covered  in  L\  (except  the  source  s)  con¬ 
tributes  to  the  weight  of  at  least  one  hyperedge  in  Li,  the  width 
W(Li)  is  no  larger  than  the  length  of  this  path  plus  one.  Also 
because  L\  is  the  shortest  path,  its  length  is  no  larger  than  the 
length  of  Lopt .  We  thus  have 

W{h)  <  C(h)  +  1  <  C(Lopt)  +  1.  (10) 

We  then  obtain  the  approximation  ratio  by  deriving  an  upper 
bound  of  C(L opt)  as  a  function  of  W(Lopt). 

Note  that  the  destination  set  Te  of  hyperedge  e*  on  Lopt 
cannot  contain  k  —  i  +  1  vertices:  its  own  source  vertex  sGi  and 
vertices  in  {se.+2,  sei+3, . . . ,  sGk,t}-  The  later  holds  because 
otherwise  Lopt  is  not  the  thinnest  path  with  minimum  number 
of  hyperedges.  We  thus  have 

k 

C(Lo pt)  <  J2(w(Lopt)  -(k-i  +  1)) 

1  =  1 

=  kW(Lopt )  -  k(k  +  l)/2 
<  W(Lopt)(W(Lopt)  -  l)/2  (11) 

where  (11)  comes  from  k  <  W(Lopt)  —  1.  Substituting  (11) 
into  (10),  we  have 

Wih)  <  W(Lopt)(W(Lopt)  -  l)/2  +  1 

<  W2(Lopt)/2  (12) 

where  (12)  holds  since  W(Lopt )  >  2. 

Based  on  (12),  if  W(Lopt )  <  y/2n,  then  W(Li)  < 

\W2{Lopt)  <  Vf  W(Lopt).  Otherwise,  we  have  W{L{)  < 
n  <  \f\A7  ( Lopt )  ■  In  summary,  SPBA  provides  a 

approximation. 

B.  For  Ring  Hyper  graphs 

Since  a  ring  hypergraph  is  a  special  directed  hypergraph,  all 
the  analysis  in  Appendix  VI- A  applies.  Specifically,  inequality 


(10)  holds.  The  problem  then  remains  in  obtaining  a  tighter 
upper  bound  of  C(Lopt)  based  on  the  geometrical  properties  of 
ring  hypergraphs. 

First,  note  that  the  length  of  a  hyperpath  L  equals  the  sum  of 
the  number  of  times  each  vertex  is  reached.  Let  Ev  denote  the 
set  of  hyperedges  on  Lopt  that  include  v  in  their  destination  sets, 
i.e., 


Ev  —  {e  £=  ^opt  :  v  E  Te}. 

Now  we  construct  a  subset  E'v  of  Ev  by  iteratively  removing 
one  from  any  pair  of  hyperedges  whose  positions  in  Lopt  are 
adjacent  until  no  such  pair  exists.  Because  at  most  half  of  the 
hyperedges  are  removed  from  Ev ,  the  size  of  E'v  is  at  least  half 
of  the  size  of  Ev ,  in  another  word  |  Ev  I  <  2\E'v\. 

Let  fimax  and  .Rmin  denote  the  largest  maximum  range  and 
the  smallest  minimum  range  among  all  vertices  in  the  given  ring 
hypergraph  Hr ,  respectively.  Let  R'min  be  the  larger  one  be¬ 
tween  Rm[n  and  the  smallest  distance  between  any  two  vertices 
in  Hr .  Based  on  the  construction  of  E'v ,  the  set  of  source  vertices 
of  hyperedges  in  E'v  satisfies  two  properties.  First,  based  on  the 
definition  of  ring  hypergraphs,  the  distance  between  any  source 
vertex  in  the  set  and  v  is  no  larger  than  the  maximum  range  of 
this  vertex  and  hence  no  larger  than  Rmax.  Second,  the  distances 
between  any  two  source  vertices  in  the  set  are  larger  than  Rmin 
and  hence  R'm[n .  Otherwise,  the  two  hyperedges  rooted  at  these 
two  vertices  can  reach  the  source  vertex  of  each  other  and  hence 
they  are  adjacent  in  Lopt  (recall  that  £  Lopt  cannot  reach  any 


vertex  in  {i 


'ei+ 2  5  sei+ 3  , 


J). 


Given  these  two  properties,  the  size  of  E’v  thus  is  upper- 
bounded  by  the  maximum  number  of  points  in  the  Euclidean 
space  that  are  at  most  fimax  away  from  v  and  at  least  R'min  apart 
from  each  other.  This  is  equivalent  to  a  sphere  packing  problem 
of  arranging  the  maximum  number  of  small  spheres  with  radius 
^min/2  inside  a  large  sphere  with  radius  Rmax  +  R'min/2.  An 
upper  bound  of  this  packing  problem  is  the  ratio  between  the 
volumes  of  the  large  and  small  spheres.  We  thus  have 

(^max  +  Rmin/^)d 


\E'\  < 


—  (1  +  2a;)0 


(^i„/2)d 

where  a  =  Rmax/ R'min-  Recall  that  \EV\  <  2\E'V\.  Note  that 
the  destination  t  can  only  be  reached  by  the  last  hyperedge 
and  hence  \Et\  =  |{e&}|  =  1.  We  thus  have 


£(Lopt)=  J2  l^l  +  l^l  (13) 

v£Lopt \{t} 

<  2(1  +  2a)d(W{Lopt)  -  1)  +  1 

<  2(1  +  2 a)dW{Lopt)  -  1.  (14) 

Substituting  (14)  into  (10).  we  have 

W{Li)  <  £(Lopt)  +  1  <  2(1  +  2 a)dW(Lopt)  (15) 

i.e.,  SPBA  provides  a  2(1  +  2a) d -approximation  for  TP  in  d-D 
ring  hypergraphs. 


C.  Asymptotic  Tightness 

We  now  prove  that  y^-ratio  is  asymptotically  tight  even  for 
2-D  disk  hypergrpahs.  The  proof  has  two  steps.  First,  we  con¬ 
struct  a  directed  hypergraph  H  for  which  the  worst  case  ratio  is 
asymptotically  reached.  Next,  we  show  a  2-D  disk  implementa¬ 
tion  of  H. 

Consider  the  the  following  hypergraph  H  illustrated  in 
Fig.  17  with  k  red  vertices  vit*..,Vk  and  k1  blue  vertices 
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Fig.  17.  Worst-case  scenario  for  SPBA. 


Fig.  18.  2-D  disk  implementation  of  the  worst-case  scenario  for  SPBA. 


ui,...,Uk'  along  with  the  source  s  and  the  destination 
t.  Each  red  vertex  Vi  has  one  outgoing  hyperedge  e  with 
Te  =  {ui, . . . ,  Vi+i}  (let  Vk+i  denote  t).  Each  blue  Ui 
has  one  outgoing  hyperedge  e  with  Te  =  {zg+i}  (let  Uk'+ 1 
denote  t ).  Finally,  we  add  two  hyperedges  that  connect  source 
s  to  vi  and  m,  respectively. 

Let  k'  =  k{k  +  l)/2  +  1.  Since  the  shortest  path  traverses 
through  the  blue  hyperedges  while  the  thinnest  path  through  the 
red  ones,  the  approximation  ratio  is  given  by 

7(*0  =  {k2  +  k  +  2)/(2k  +  4).  (16) 

Note  that  the  total  number  of  vertices  is 

n  =  k-\-k'-\-2  =  k  +  k(k  +  l)/2  +  1. 

When  n  is  large,  k  ~  \f2n  and  a(k)  ~  |  ~ 

Next,  we  implement  the  above  hypergraph  under  a  2-D  disk 
model  as  illustrated  in  Fig.  18.  The  red  vertices  are  located  on 
a  straight  line  with  RVl  =  RV2  =  1,  Rv.  =  2l~2  for  i  >  2. 
The  source  vertex  s  is  located  on  the  line  to  the  left  of  zq,  and 
both  its  maximum  range  and  its  distance  to  v\  equal  RVk .  The 
terminal  vertex  t  has  a  maximum  range  of  0  and  is  located  to 
the  right  of  Vk  with  a  distance  of  RVk .  The  maximum  range  of 
a  blue  vertex  zq  is  RVk  —  ie  where  e  is  a  small  positive  value 
to  prevent  zq_i  from  overhearing  messages  transmitted  by  zq. 
The  blue  vertices  are  located  on  a  route  from  s  to  t  that  contains 
two  vertical  line  segments  of  length  (1  +  l)RVk  and  a  horizontal 
one  of  length  3 RVk ,  as  demonstrated  by  the  blue  dashed  lines 
in  Fig.  18.  The  positive  parameter  l  is  used  to  prevent  a  blue 
vertex  from  overhearing  the  last  red  vertex  Vk .  In  the  asymptotic 
regime  with  large  k ,  l  can  be  set  sufficiently  large  so  that  the  k' 
blue  vertices  can  be  implemented  along  the  depicted  route  from 
s  to  t. 

Appendix  G 
Proof  of  Theorem  8 

Let  L2  denote  the  path  in  hypergraph  H  from  s  to  t  given  by 
the  TSBA  algorithm  and  Lopt  the  thinnest  path.  Let  Li(v)  and 
L2  (v)  denote  the  paths  from  s  to  a  vertex  v  given  by  SPBA  and 
TSBA,  respectively.  The  following  lemma  establishes  a  prop¬ 
erty  of  L2(v). 

Lemma  7:  For  any  hyperedge  e  in  H,  we  have,  \/v  G  Te 
W(L2(v))  <  \L2(se)UTe\. 

Proof:  Lemma  7  follows  directly  from  the  tree  structure  of 
TSBA.  ■ 


A.  For  General  Directed  Hypergraphs 

Let  L0 pt  =  {ei, . . . ,  e^}  denote  the  thinnest  path.  For  ease 
of  presentation,  let  the  sequence  of  source  vertices  sei , . . . ,  sek 
and  the  final  destination  t  be  denoted  as  vi, . . . ,  Vk+i-  Let 
U  =  •  Based  on  Lemma  7,  we  have,  for  allz  =  1, . . . ,  k 

W{L2(vi+1))  <  \L2(vi)UTe.\ 

<  W (. L2(vi ))  +  |Tei\{ui, . . . ,  Ui+i}|  +  1 

=  W(L2(vi))  +  \Tei\U\  +  l  (17) 

where  (17)  holds  since  Tei  does  not  contain  vertices  in 
{sei+2, . . . ,  sefc,  t}.  Summing  (17)  over  z,  and  noticing  that 
W(L2(v i))  =  1  and  L2(u^+i)  =  L2,we  have 

k 

W(L2)<k  +  l  +  J2\TeAU\. 

i— 1 

Next,  since 

W(Lopt)  =  k+  1  +  juti  (TeiW)| 
we  can  upper-bound  \Te.\U\  by  W (Lopt)  —  k  —  1  for  any  i  = 
1, . . . ,  k.  Thus 

W{L2)  <  k  +  1  +  k(W(Lopt)  -k-  1). 

The  right  side  of  this  inequality  is  a  quadratic  function  of  k  with 
the  maximum  at  k  =  W (Lopt)/2.  We  thus  have 

W(L2)  <  1  +  W2{Lopt)/ 4. 

If  W (Lopt)  <  2 \fn  —  1,  the  approximation  ratio  7  is  given 
by 

7  =  W(Lopt)/ 4  +  l/W(Lopt)  <  n/ 2VT^T.  (18) 

The  inequality  holds  because  the  function  |  +  Ms  an  increasing 
function  for  x  >  2. 

IffL(Lopt)  >  2 y/n—  1,  we  have 

7  <  n/W{Lopt)  <  n/{2V^l).  (19) 

This  completes  the  proof  for  case  of  general  directed  hyper¬ 
graphs. 


B.  For  Ring  Hypergraphs 


Let  L\ (■ t )  =  {ei,  e2, . . . ,  e^}  be  the  shortest  path  from  s  to  t. 
Let  Vi  denote  the  source  vertex  of  e*  and  u&+i  =  t.  We  prove, 
through  induction,  the  following  inequality  for  all  i  =  1 , . . . , 
k  +  1: 


W(L2(vl))<£(L1(vi))  + 1.  (20) 


When  i  =  1,  (20)  holds  since 


W{L2(v1))  =  1  C{L1{v1))  =  0. 

Now  assume  that  (20)  holds  for  i  —  1,  i.e.,  W(L2(vi- 1))  < 
£(Li(vi-i))  +  l.  Based  on  Lemma  7  and  this  induction  assump¬ 
tion,  we  have 


W(L2(Vi))  < 


^2(^-1)  U  Tei_1 


<W{L2{vi^))F\Tei_^ 

<  C(Li(vi-i))  +  1  +  | Te._1 1 

=  C(L1(vi))F  1. 

This  completes  the  induction.  Considering  z^+i  =  t,  we  have 
W{L2)  <£(Li(t))  +  l.  (21) 

From  (15),  (10),  and  (21),  we  have  W(L2(t))  <  2(1  + 
2a)dW(Lopt),  i.e.,  TSBA  provides  a  2(1  +  2a)2  approxima¬ 
tion  for  ring  hypergrpahs. 


C.  Asymptotic  Tightness 

We  first  construct  a  directed  hypergraph  H  as  illustrated  in 
Fig.  19.  The  vertex  set  of  H  consists  of  two  types  of  vertices: 
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k  +  1  normal  vertices  vq,  . . .  ,Vk~i  and  t,  and  k  super  vertices 
u\ , . . . ,  Uk- 1 ,  u,  each  containing  &  —  1  normal  vertices.  Rooted 
at  each  normal  vertex  Vi  are  two  hyperedges  e*+i  and  e'+1 .  Hy¬ 
peredge  ei  has  destination  vertices  Te.  =  {vi,Ui}  and  hyper¬ 
edge  e\  has  destination  vertices  Te>  =  {17,  u}. 

It  is  easy  to  see  that  the  thinnest  path  from  vo  to  t  is  Lopt  = 
{e^ ,  e'2 , . . . ,  e'k ,  e'k}  with  width  2k.  However,  TSBA  returns  the 
path  Ld  =  {ei ,  e2 , . . . ,  e/t }  with  width  k2  + 1  in  the  worst  case.9 
The  approximation  ratio  is 

W(L)  _  k2  +  1  _  n 
7  “  W(Lop t)  ~~  2fc  ~~ 

Given  the  similarity  between  H  and  the  hypergraph  Hi  con¬ 
structed  in  the  proof  of  Theorem  1 ,  we  can  follow  the  same 
approach  given  in  Appendix  A  to  implement  H  under  a  2-D 
disk  model.  However,  this  implementation  requires  additional 
vertices  (referred  to  as  auxiliary  vertices)  that  may  render  our 
previous  approximation  analysis  invalid.  To  maintain  the  ratio, 
each  original  vertex  (including  the  vertices  in  a  super  vertex)  in 
H  is  replaced  with  c  vertices  (clustered  together)  in  its  2-D  disk 
implementation,  where  c  is  the  number  of  auxiliary  vertices  in¬ 
troduced  by  the  implementation.  In  this  case,  TSBA  returns  a 
path  that  covers  {ui, . . . ,  Uk-i,  vq,  . . . ,  Vk}  along  with  a  set  of 
auxiliary  vertices.  The  thinnest  path  covers  {u,  uo, . . . ,  Vk}  and 
another  set  of  auxiliary  vertices.  The  approximation  ratio  in  this 
2-D  disk  hypergraph  is  given  by 

_  (fc2  +  l)c  +  c'  _  fc2  +  1  +  j- 
7  2  kc  +  c"  2  k  +  — 

C 

where  c'  and  c"  denote  the  number  of  auxiliary  vertices  covered 
by  the  path  returned  by  TSBA  and  the  thinnest  path.  Since  ^  < 
1  and  —  <  1,  when  n  is  large,  we  have  7  ~  0  A— r,  i.e.,  the 
approximation  ratio  2^— _  is  asymptotically  tight  in  2-D  disk 
hypergraphs. 
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